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DEVELOPMENT OF DATA PROCESSING, 
INTERPRETATION AND ANALYSIS SYSTEM 
FOR THE REMOTE SENSING OF TRACE 
ATMOSPHERIC GAS SPECIES 


INTRODUCTION 

The Cooperative Agreement, NCC1-34-, between NASA/Langley 
Research Center and Old Dominion University, entitled "Develop- 
ment of Data Processing - Interpretation and Analysis System 
for Remote Sensing of Trace Atmospheric Gas Species" has 
represented a major research effort during the period of 
January 19'81 to November 19 81, Research work performed by ODU 
specifically addressed problems associated with the development 
of the MAPS experiment program being conducted at NASA/LaRC. 

The primary thrust of this research has been the utilization 
of ’the MAPS experiment data in three application areas: 

(1) low altitude aircraft flights (one to six km) 

(2) mid-altitude aircraft flights (eight to 12 km) 

(3) orbiting space platforms. 

Extensive research work iu four major areas of data 
management has been the framework for implementation of the 
MAPS experiment technique. These areas- are: 

(1) data acquisition 

(2) data processing, analysis and interpretation 
algorithms 

(3) data display techniques, and 
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(4) information production. 

The accomplishment of these tasks under this agreement 
Is ongoing and has required major developmental research and 
analysis of two la-rge- MAPS experiment data sets, the 1978 St.. 
Petersburg-, Florida data and the. 19 79 MONEX data. Portions 
of the methodology and operational software developed, in the 
analysis of these two extensive- sets of data have been applied 
in -the MAPS/OS, TA Space Shuttle data reduction plan and soft- 
ware schema. 

A major contribu.ting factor in the successful progress 
of these research tasks has been the implementation of -the 
Tektronix 4052. desktop graphics system in -the computation work 
requirement.. Use of the Tektronix 405'2. desktop graphics, 
system, bo-fch in the stand-alone mode and as, an online, terminal, 
increased significantly during this 'reporting period. This 
system has shortened the turn-around time for computer results., 
especially in the form of graphics output, and has added to our 
capability of summarizing data with quickly produced,, highly 
accura.te graphs and with mathematical curve fits . Diverting 
computer work, to the Tektronix system has been a contributing 
factor in the reduction of expenditures, approximately a 
factor of two savings' by the MAPS ODU personnel, for "the use 
of NASA operated computers. Currently, computer costs are 
wi-t±iin budget rather than exceeding allocations as has occurred 
in recent fiscal years. 

Large amounts of data are currently stored on Tektronix 
floppy disks at essentially no cost, discounting -the initial 
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investment in this computer system. If required/ this data, 
can be transferred to the NASA computing facility, a process 
which is not excessively time-consuming when accomplished at 
1200 baud. This equipment is located off-site in. the Research 
Drive' office complex and is. accessible at all times, therefore 
grea-tly extending our- capabilities beyond the limlta.tions of 
regular weekday accessibility to NASA computing facilities- 
"Down time" due to hardware failure has been negligible, and 
hardware data processing has been excellent,. 

Operating as an on-line hookup to the NASA compater net- 
work, -the Tek,tronix system has grea.tly decreased -the time 
required to "debug" new graphics software.. Past unavoidable 
delays in -the generation and • delivery of paper plo.ts were, 
eliminated when graphics output could be previewed, on the 405.2 
system terminal' screen. Corrections., changes, and enhancements 
to the graphics images' are now possible before the final, output, 
is routed to a permanent medium plotting device. High quality 
graphics output are produced by the Tektronix 4663 interactive 
digital plotter. These images are easily enlarged or reduced 
according to -the needs of -the' user.. Multicolor plots can also 
be generated to add emphasis, clarity, or simple eye-appeal. 

The DEC- 10 computer operated at the main campus of ODU is also 
accessible to the Tektronix system. To facilitate program 
development, files such as the meteorological models under 
development were transferred from ODU to NASA and vice versa 
via telephone lines wi-th intermediate storage on -the Tektronix 
cartridge tape. 
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Offline this system was used largely for data display and 
curve fitting, using both software purchased from the vendor 
and code created in-house. The BASIC programming language 
built into the system was mastered readily by all users. In 
addition to coding numerous special purpose programs in BASIC, 
a general utility plotting package was written and executed 
extensively. This' program displays one or more sets of 
coordinate pairs- on either linear or logarithmic axes.. The 
user controls' scale factors, titles-, symbols-, and -the- optional 
inclusion of a curve representing a "best fit" of -the data by 
one of several linear' equa,tions. This package also displays 
statistics' including maximum and minimum residuals and esti- 
mated values- of the dependent variable. The implementation of 
this- computer graphics system in the performance of research 
work by ODU personnel has' resulted in Incalculable, cost savings' 
in computer expenditures by NASA. 

Although the research being performed by ODU is currently 
under a continuation agreement of NCCl-34, specific develop- 
mental accomplishments in the data management organization have 
been successfully performed for low altitude and mid-altitude 
aircraft experiments. These research tasks addressed work in 
three areas, St. Petersburg Test Flight data, the 1979 Summer 
MONEX data, and MAPS/OSTA Experiment integration and testing. 
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ST'.. P.ETERSBURG FLIGHT. TEST - LOW ALTITUDE 

In August 19 7-8', low alti.tude aircraft flight tests were 
conducted using' the MAPS Brassboard GFCR to measure the 
column burden of carbon monoxide (CO) in the St. Petersburg - 
Cleaor'water urban areas of Florida. The data acquired during 
these flight tests provided a data- base- for the evalua.tion of. 
the performance of the GFCR sensor system' and for the testing 
of refinements to the data, reduction me-thod of the MAES experi,— 
ment. Research work, performed utilizing this GFCR data improved' 
the data reduc.tion. procedure for MAPS.. 

These, flight tests;- evaluated the effects of the MAES 

instrument output signal to hot-and humid sub.tropical a.-tmos- 

pheric- conditions;' to partial cloud cover in the- field of view-;. 

and. to' varying, surface- (land or: water) emissivities-. Analysis' 

\ 

of independently obtained CO vertical mixing ratio profile data- • 
obtained from -fche flight tests was- used to characterize the 
spatial and temporal variations of CO within this urban study 
area. 

The St. Petersburg Peninsula on the west coast of Florida 
was selected as a. test site' because it is loca.ted in a sub~- 
tropical climate and is adjacent to the Gulf of Mexico, which 
is a non-source CO area. This test site also provided a 
data se.t with variable cloud distribution (both type and 
coverage) as well as variable contamination sources. Sufficient 
emission inventory and simple geographic distribution of sources 
at this test site simplified the data analysis procedure used 
in the verification of the MAPS CO measurement technique . 
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Experimental Details 

To perform the St. Petersburg' experiment, a Cessna 40 2B 
(N69 402) twin engine, aircraft wi-th a crew of -three {pilot, 
flight coordinator, and sensor, operator) was. equipped wi.-th 
•the MAPS Brassboard GFCR, Precision Radiometric Thermometer 
(PRT-5) , hygrometer system,, total temperature probe, altitude 
pressure transducer, Loran route verification system, 3-5' mm 
aerial camera., and a Loran navigational, sys.tem. 

Six flights were flown from August L4' to August 1.7, 19 78, 
over -the St. Petersburg - Clearwater area, as shown in Figure 
L. Remote sensing data were collected during, three flights 
at an altitude of 6.1 km and three flights at an altitude of 
2. 4' km. Shown in Table 1 are the beginning and ending' times 
of the flights and the latitudes at which the data, were 
collected. The first and third flights' were chosen for. 
de.tailed data reduction and, will, be referred to in -this paper 
as Flight I and. Flight II, respectively. 

For each flight the Cessna 402B departed the St, Petersburg 
Clearwater International Airport, ascending 'at a rate- of 3.6 
m/sec to flight altitude. This procedure required about 30. 
minutes during . which time the MAPS instrument was cycled 
through a balance/calibration sequence. The aircraft pilot 
entered -fche test pattern at location BALL labeled in Figure 1 
and proceeded in a counterclockwise flight direction along 
the defined aircraft ground track, requiring approximately 45 
minutes . 
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Figure 1. St. Petersburg ground track. 
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Table T 

Beginning and ending flight times 


Flight no. Date 


Approximate 

Start Stop 


Alt 


km 


1 

2 

3 

4. 

5' 


Aug. 14, 1978 1458 hours (LOT) 1740 hours.(LDT) 


Aug. 14-, 1978 1828 2.100 

.Aug. 14, 1978 2152. 0100 

Aug. 1‘5, 1978 ' 0125' 0430 

Aug,. 16, 1978 0614- 0930 

Aug. 17, 1978 1018 


6.1 


2.4 

6 . 1 ' 


2.4 ' 


6.1 


6 


1253, 


2.4 
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The MAPS instrument was recycled through a balance/cali- 
bration sequence prior to the aircraft commencing a race track 
pattern, descending over BALL. During the descent, the air- 
craft was flown straight and level for approximately two 
minutes at each of the nine different altitudes where GFCR 
data and ambient air samples were obtained. 


Measurements’ 

The inflight data collected consisted of radiometric 
surface temperatures; two GFCR output signals: ^^cell^ 

containing the gas of interest, and a vacuum cell.; 

ambient air samples, total air temperature at flight altitude, 
dewpoint temperatures at flight altitude, and altitude in 
feet. 


Inflight, analog output signals of -fche flight instrumenta- 
tion was recorded by a Pulse Coded Modulation (PCM) digital 
recording system. The analog signals of the instruments were 
digitized and stored in binary format on seven- track computer 
compatible tapes. (1) The PCM sys.tem recorded one record 
every twelve seconds. Records contained 31 frames of data and 
were separated by a one frame inner record gap . 


Data Reduction Procedure 

A schematic of the procedure used to reduce the GFCR data 
obtained at St. Petersburg is shown in Figure 2. The first 
step in the data reduction procedure was to utilize computer 
code to select every fifth frame of data from the PCM tape. 



original page is 



Figure 2. GFCR data reduction procedure 
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starting with the firs^t frame, and create- a file- which 
retained the same 32-frames per record format as the original 
tape. This file, which contained only a sampling of the PCM 
data, was used as input to another computer program which 
averaged selected parameters within each frame. The averaged 
data corresponds to a • frame sampling rate of approximately 
2.67 frames, per second. A flow diagram of this procedure is 
shown, in Figure- 3 . 

The data processing procedure to this' stage was 
accomplished on the CDC Cyber computer network maintained by 
the Analysis and Computation Division (ACD.) of NASA/L'aRC. 
Averaged values of •the data were transferred to the Tektronics 
4052 graphics system at -the Old Dominion University Research. 
Foundation office i,n Hampton for' final, data, reduction. 

The ‘next step in the data reduction process was to use. 
the meteorological parameters, outside air temperature, 
pressure, and dewpoint temperature,' obtained during -the 

t 

periodic descent of the aircraft, .in the initialization of a 
computer program that simulated- -the overall a-bmospheric 
radiation characteristics during the time of -the GFCR measure- 
ments. This modeled atmosphere, which characterized the 
temperature, pressure., and water .vapor distribution, was used 
to compute the optical -thickness specifications of a multi- 
layered homogeneous a-fcmospheric model. This" modeled atmosphere, 
composed of -the .average layer values for the temperature, 
pressure, and volumetric water vapor mixing ratios, was used 
as the initial atmospheric conditions for SMART.' A major 
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Figure 3. PCM data reduction procedure 
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element in the data reduction procedure/ the radiative 
transfer calculations, was accomplished by using the SMART 
program. The tables of theoretical radiances (N) and 
differential radiances (AN-) generated by SMART ware 
parameterized, thus establishing a mechanism by which CO 
mixing ratios could be obtained for the selected atmospheric 
conditions. The next step was to use the GFCR flight data 
and the parameterized SMART values as input into a computer 
program that computed the CO mixing ratio for a given radiance 
(N) and differential radiance CAN) value. Plotting CO mixing 
ratios versus time in flight was the final step. 

One input parameter to the SMART program is the number 
of homogeneous absorbing layers, which comprises the total 
number of- modeled atmospheric, instrument, and calibration 
layers. The modeled atmosphere was partitioned at the 
descent altitudes of the race track pattern, ensuring no large 
variation in the modeled atmospheric temperatures within the 
vertical dimension of a layer. 

The CO mixing ratio values were chosen to bracket -the 
expected range of CO values in -fche geographical area of 
interest. Shown in Table 2 are -fche mixing ratios obtained 
from the gas multipliers used for this study. The background 
CO mixing ratio was assumed to be uniformly distributed in the 
vertical wi-th a CO mixing ratio of l.OE-7. 

Shown in Table 3 are the specified model atmospheres for 
Flight I and Flight II. The temperatures of the layers are 
expressed in degrees Kelvin, the pressure of each layer is 
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Table 2 

Mixing ratio values and surface temperature 


CO mixing ratios (ppbv) 


0 

150 

50 

200 

75 

300 

100 

400 

125 

500 


Surface temperatures (K) 


298.16 

308.16 

300.16 

310.16 

302.16 

312.16 

304.16 

314.16 


306.16 


316.16 
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Table. 3 

The; model ed atmospheres for Flights I and IT 


Temperature- 

K 


Pressure 
atm ' 


Thickness 

cm 



Flight I 


302.6 

.9890 

15200. 

301.9 

. 9720. 

16800. 

299.3 

.9440 

32000. 

296.9 

.,9090' 

32000. 

294.1 

.8760 

29000. 

292.3 

.8450 

30400'. 

290 .T 

.8,150 

29000. 

288.4- 

.7850 

32000. 

285.5 

.7550. 

30500. 

283.7 

.7140 

57900. 

279.6 

.6370 

121900. 

274-. 3 

.5650 

6.1000. 

269.9 

.5210 

60900, 

266 . a 

.4890 

30500. 

266-..0 

.4680' 

33500.. 


Flight IT 


300.50 

.9890 

15200, 

299. 60- 

-.9720- 

16800. 

298 . 20 

.9440 

32000. 

295-. 90 

.9090- 

32000. 

292.90 

.8760 

29000, 

291.40 

.8450 

30400. 

290.60 

.8150 

29000. 

288.90 

.7850. 

32000. 

286.60 

.7550 

30500. 

285.10 

.7140 

57900. 

280.60 

.6370 

121900. 

274.90 

.5650 

61000. 

279.20 

.5210 

60900. 

266.20 

.4890 

30500. 

266.-40 

.4680 

33500. 
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expressed in atmospheres/' .ajad. the thickness of each layer is 
expressed in centimeters - 

The parameterization of the theoretical radiance and 
dif fe-rential radiance values computed in SMART was accomplished 
using mathematical fits to represent the- CO mixing ratio as a 
bi-variate function (i.,e. , the CO mixing ratios are a function 
of radiance and differential radiance values) . The surface 
temperatures spanned tile range of measured surface tempera- 
tures, observed during the flights,, 

Mathemati-cal fits were- used to parameterize the bivariate, 
characteristics of the SMART radiance values. A fourth degree 
least squares polynomial fit was used to represent CO mixing 
ratios- as a function of.’ theoretical differential radiance (AN) 
values for a famiLy of varying surface temperature curves. 

A third order least squares polynomial fit was used to represen.t 
the coefficients (a) obtained from the first fit as a function 
of a particular theoretical radiance (N) value. The radiance 
value used for this study was based on a nominal CO mixing 
ratio of 100 ppbv. The coefficients (3*s) obtained from the 
fit of the first coefficient will be used to determine the 
mixing ratios of CO from the GFCR data. . 

The inflight GFCR data were converted from electrical 
units to physical units by linearly interpolating on the cali- 
bration curves'. The radiometric and differential radiance 
signals in the bivariate fits were used to determine the CO 
mixing ratios by interpolation. This was accomplished by 
first determining the a*s from the radiance values and the B 
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coefficients'. Then the a * s^ and differential radiance are. 
used in an equation -that infers the CO mixing ratios for a. 
particular radiance and differential radiance value.. The 
computations sho'wn below use a. radiance (N} and differential 
radiance (DN) obtained from inflight data. The beta 
■coefficients were computed by parameterization of SMART 
tables generated for this- flight. The. alpha coefficients. 
(Eguatipn. set II.) were- computed from the. given B and. N' values, 
and used, along with DN, to determine the mixing ratios-. Data, 
were, taken at an altitude of 6-..1 km over a water surface. 

Relative, differences between fit. and corresponding 
theoretical SMART values are summarized- in, Table 4. These, 
data are applicable to wa.ter and land surfaces at temperatures 
of 29'8°K and 304-^K. Figure 4- is a plot of the- bivariate fit 
(solid curves) and the SMART discrete values (symbols) 
representing da.ta values over a, water, surface at three 
surface temperatures (298°K, 304°K, and 306°K) . 
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MR-, = a-,., + a .,2 ON + a -,3 DN^ + DN^ 

MR2 = a 2 i “22 ‘^23 “24 

MR^ = a2i +• a^2 ON +- DN^ + DN^ 

^^4 =■ “4.1 ■^' “42 “43 “44 DN^ 

*^^5 '*■ “52 “53 ‘^54. 

^^6 =• “61 '■' “62 “63 ^ “64 


“xl ^11 ^12 ^13 

“x2 ^ ^12 ■^' ^22; ^ ^23 

“x3 ^ ^31 ■'■■ ^32 ^.33.^^ 

“x4 ^41 ^ ^42 ^43 


( 1 ) 


(ID 
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N = 2.52E-5 (watts/cm^ - sr) 
DN = 4.T5E-7 (watts/cm^ - sr) 


4.627E-5 
1.838E-10 
6.205E-14 

= -2.09E-8 
= .397 

=■ -5.32 X 10^ 

=• 1 ..04 


Beta Coe.ff1cients 
9.832E.-8 7. 733 Erl 

-6.419E-3 -2.017E-4 

6.694E.-1 2.077E.-8 


9.441E-12 
-5.202E-.17 
7.4-1 ie;-21 


Mixing ratio = 126.4- ppbv 
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Table^ 4. 

Errors associated with mathematical fits for water and 
land surfaces at 6.1 km 


Theoretical 
mixing ratios 
ppbv 


Fit 

mixing ratios 
ppbv 


Relative 

difference 

% 


50 

48 

-2.6 

8 = .88 

75 

100 

73 

99 

-2,.0 

-4.4- 

Tsurf = 2S8 

125 

T24 

-.6 


150 

151 

0.4 


200 

205 

2.5, 



50 

47 

-5.8 

e = .88 

75 : 

100 

72 

98 

-3.4 
-2,. 5 


125 

124 

1.0 


150 

151 

0.8 


200 

210 

5.0 



50 

48 

-3.1 

e =■ .98 

75 

100 

74 

99 

-1.6 

-1.0 

l“surf = 298 

125 

125 

-.3 


150 

151 

.7 


200 

206 

3.1 



50 

47 

-5.2 

£ = .98 

75 

73 

-2.7 

’’’surf " 

100 

98 

-2.0 

125 

125 

-.6 


150 

150 

1.3 


200 

200 

5.8 






CO MIXING RATIO ♦ lOE-7 


Figure 4. plot of fit values versus discrete values 
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DATA RESULTS AND AISIALYSIS 

To completely analyze the GFCR data obtained over St. 
Petersburg/ it was necessary to characterize the prevailing 
meteorological conditions during the period of August 13 
through 18, 1978, in the urban area of St. Petersburg, 
Meteorological data were obtained from several sources.* The 
850 mb and 500 mb analyses were supplied by the National 
Weather Service. (NWS) . Specially launched rawinsondes_, 
corresponding to the approximate times of the aircraft 
descents, were obtained from the NOAA station at Ruskin, 
Florida. Normal rawinsonde data for Ruskin and surface 
observations for Tampa International Airport, St. Petersburg — 
Clearwater International Airport, Albert Whiting Airport, and 
MacDill Air Force- Base were obtained from .the National 
Climatic Center (NCC) in Asheville, North Carolina. The out- 
side air temperatures and dewpoint temperatures were obtained 
from aircraft data. Satellite imagery was provided by NOAA. 

The atmospheric profiles obtained from the aircraft 
(Figure 5) were compared to the atmospheric profiles obtained 
from the rawinsonde data (Figure 6) for Ruskin, and minor 
differences were observed in the temperature profiles for both 
flights studied. The implications associated with the n\inor 
differences have 'minimal effect upon the interpretation of the 
atmospheric profiles. These differences in the meteorological 
profiles can be attributed to differing sites and times of 
data acquisition. Aircraft data were taken over water, and 
rawinsonde data were acquired in land. These comparisons 
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Figure 5* Aircraft profiles for Flights I and II. 
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FLIGHT I 



TEWPERATURE (C) 



FLIGHT II 



TEMPCRATURE (C) 

Figure -Rawinsonde versus aircraft temperature profiles. 
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validate the aircraft data and increase the feasibility of 
using the upper, level winds reported by Ruskin in an inte- 
grated, detailed analysis of the meteorological, and CO data.. 

CO Measurements 

The first data collecting mission was flown on the l-4th, 
commencing at approximately 3:00 P.M. (LOT). The upper level 
winds at 6,1 km and 5 ..5 km were out of the south- southwest 
between 2 and 4-..6' m/sec,, respectively. Figure 7 shows the 
GPCR CO mixing ratio measurements and upper level W'inds 
superimposed on a map of the area... The average measured CO 

mixing ratio on the easternmost data leg was 142 ppbv, on the. 

* 

northernmost data leg' was 135 ppbv, on, the westernmos.t data 
leg Cover water) was 112 ppbv, and on the southernmost data 
leg was. 112 ppbv. 

The second flight started at about 9:00 P.M. on the 14th 
with the winds out of the north-northwest at approximately 
1 m/sec. Shown in Figure 8 are the CO mixing ratios and 
upper level winds sketched on a background of the region. 

The average measured CO mixing ratio on the eastern data leg 
was 86 ppbv, on the northern data leg was 83 ppbv, on the 
western data leg was 89 ppbv, and on the southern data leg 
was 89 ppbv. 

Vertical profiles of CO were obtained over BiALL for both 
flights after the constant altitude data circuits had been 
completed. Table 5 contains CO values measured by a gas 
chromatograph system *and CO values from the GF.CR analyzed data- 
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Figure 8- CO mixing ratios and upper level winds for Fliohf- tt 





Vertical CO mixing 




Altitude 

km 

GC Flight I 
mixing ratio 
ppbv 

- 

Point 

Col umn 
burden 


r 

6.1 

90 

107' 

5.5 

80 

no 

4.3 

129 

117 

3.0 

115 

114 

2.4 

- 

- 

1.8 

95 

113 

1.2 

105 

125 

0.6 

143 

157 

0,3 

- 

- 

0.1 

226 

226 
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ble 5. 

ratios for 6C and GFCR 


GC Flight II GFCR Flight II 
mixing ratio mixing ratio 


ppbv 


Point 

Col umn 
burden 

70 

84- 

73 

86 

69 

88 

77 

94 

112 

100 

88 

94 

90 

99 

106 

' 115 

140 

140 


ppbv 
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The independent gas chromatograph system has been used in 
previous CO- measurement analysis and has been proven reliable. 

(2, 3, 4). The gas chromatographic data represent an in situ 
point source measurement while the. GFCR represents a, column 
burden measurement. To relate- the data obtained by’ the 
in situ analysis and the- GFCR, it was necessary to calculate 
an integrated atmospheric, column burden bounded, by the height 
of the. sensor and the surface., using multiple ambient, air 
sample bottle gas' chromatographic- measurements-. 

The GFCR and GC mixing ratios were compared at flight 
altitude and at -several levels during the descent cf. the air- 
craft for the second flight. The comparison at. flight altitude 
(6,. l. hm) was a. column burden measurement comparison, and a. 
relative difference of 4-%‘ was- compu.ted,. K' secondary comparison, 
made- at several des.cent levels is- shown graphically in Figure 9 .• 
The relative differences in measurements range from 2.% to 22.%. 

For the- first flight the two systems were only compared 
at a flight altitude of 6.1 km. The results of this column 
burden measurement comparison indicated a relative difference, 
of approximately 18%. 

The large relative differences in -the CO data measurements 
obtained during the descent of the second flight and at flight 
altitude 'of Flight I are. a result of the same general factors. 

In bo-th cases the radiative' transfer calculations were per- 
formed using a CO profile uniformly distributed. However, the 
analysis of the GC CO data as shown -in Figure 9 clearly shows 
that this condition was nonexistent in the atmosphere at the. 
time- of the CO measurements. In particular, gradients of CO 
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were shown in the middle of the air. column. An in-depth 
analysis of. the temperature profile for the descent shows a 
temperature inversion to be located at approximately 2.5 km. 
Acting as a cap, this temperature inversion produced the CO 
gradient in the middle portion of the air column, which could 
result in a decrease in the accuracy of the measurements 
obtained by the MAPS GFCR. Therefore, the results' of the 
GFCR and GC CO data comparison show that the GFCR CO values 
for these two flights are very reasonable.. 

Box Model 

The CO data obtained from these two independent systems 
and, the prevailing- meteorological conditions can. be used in . 
con.junction with a CO emission inventory in applying a. con- 
servation of mass model to the. study area. In order to 
perform these conservation of mass calculations', a source 
estimate of CO is required. This estimate was obtained from 
16 traffic counting stations which were strategically set up 
to monitor the flow of vehicular traffic into and out of -the 
study area.- The station locations are shown in Figure 10. 

To perform the mass conservation calculations utilizing 
the data described above the Flight I in the St. Petersburg 
area, the following assumptions are made: 

- Sources are uniformly distributed. (Determined from 
the highway network. ) 

- Rate of emission of pollutant per unit area is constant. 

- The net diffusion of pollutants to the sides are 
neglected. (The winds did not vary much oyer the time 




Figure 10. Traffic counter locations. 
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of interest. 

- A series of closely spaced line sources was assumed. 
(Determined from highway network.) 

— Only automobile emissions are to be considered. (No 
industrial sources .inside area of interest..) 

- Automobile emission of' 112. g/day was used.. (5) 

Between Flight I and Flight II the a.tmosphe-ric conditions' 

transitioned such -that for -Flight II several of the assimptions 
required to perform -the mass conservation calculations for an 
area source were violated. Consequently, the results of the 
mass conservation calculations for the second flight were not 
performed.. 

To compute the expected CO mixing ratio for the first 
flight, given the previousiy mentioned assumptions, Equa.tion 
(1) was used. (6) . 


X 


2 X (Q/L^) 

(7T . u . D )^-^^ 
2 


^ 1/2 


Cl)' 


Where x the computed concentration, Q is -the emission rate, 
L is the leng-fch of a side of a surface area of a square, u is 
the average wind speed, X is -hhe- distance to be traversed, 
and is the eddy diffusivity. 

The volumetric mixing ratio was computed according to 
Equation (2) 

Q„ ^ 22.4 X 10“^ (Q /M ) . (T/273) (2) 

V m p 

where Q and M are the mass emissions and molecular weight, 
m p 

respectively, of the gas species of interest, and T is the 


surface temperature. 
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The results of the box model calculations for Flight I 
show an increase in the mixing ratios of CO of approximately 
35 ppbv as the source area is traversed from point A to. point 
B (Figure 11) and a 28' ppbv increase from point C to point D 
(Figure 11) . The computed increase is approximately 22% 
higher than the background mixing ratios of CO measured by 
the GPCR' over the eastern portion of the ground track', at an. 
altitude- of approximately 6.1 kiti. 





Figure 11. Area of the box model. 
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SUMMER MONEX - MID-ALTITUDE 

The-MONEX (MONsoon Experiment) of 1979 represented one 
of the most intensive scientific measurement projects which 
investigated the a.tmospheric interhemi spheric transport in 
the Indian Ocean area. The research performed utilizing the 
MAPS CO experiment measurements and the GC measurements, 
represents a very complex set of data obtained from a mid- 
altitude flight program. Specific tasks accomplished 
include: 

- production of a comprehensive set of meteorological 
data to Dr. Reginald Newell of M.I.T. and Ms. Estelle Condon 
of NASA/Ames Research Center. This data was used to analyze 
the atmospheric interhemispheric transport of carbon monoxide 
during the MONEX test period. A publication on the results 
of this analysis is currently in print. 

- categorization and sorting of all MAPS radiometer data 
for efficient data reduction. This task included the develop- 
ment and performan'ce of background target sorting schema. 

t 

In addition, calibration analysis was automated, decreasing, 
by a factor of four, the time required for processing large 
numbers of calibration data sets. 

- verification of meteorological models utilized for the 
microprocessor real time data reduction. 

- development of a more efficient and faster comprehensive 
data reduction methodology for the GFCR post-data reduction. 
This resulted in an approximate factor of 75 reduction in the 
cost of the interpolation process . 
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MAPS/OSTA-1 EXPERIMENT 

Research tasks during this reporting period have 
included major organizational, developmental, and operational 
work in support of the MAPS/OSTA Experiment successfully 
flown on STS-2 in November 1981. This work involved con- 
siderable time and effort in pre-mission and mission planning 
and operations. The major task areas of work were: 

- performance of STS—2/0STA— 1 MAPS Experiment signal 
interpretation, calibration and integration testing. 

- providing meteorological support data for aircraft 
and MAPS Space Shuttle associated operations, 

- developing a preliminary data management plan for 
MAPS/OSTA data reduction.' 

Comprehensive information on the details of the work 
performed in these- task areas will be forthcoming in technical 
reports on each of the specific accomplishments. 
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